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I . INTRODUCTION 


Our  previous  years  of  continuous  effort  (as  summarized  in 
Reference  1)  have  resulted  in  significant  progress  in  analyzing 
the  flow  field  over  powered  missiles,  especially  at  transonic 
speeds.  It  is  important,  in  addition  to  our  previous  efforts, 
to  estimate  influence  by  fins,  surface  cutouts  and  different 
incidence  angle  ranges,  etc.,  to  the  missile  performance.  In 
order  to  obtain  an  inviscid/viscous  flow  interaction  solution, 
a complete  understanding  of  the  potential  flow  field  about  mis- 
sile configuration  is  required.  We  have  worked  toward  this 
goal.  Nevertheless,  the  originally  planned  three  years'  effort 
was  terminated  at  the  second  year  due  to  the  change  in  funding 
situation.  During  the  period  work  performed,  three  papers  were 
prepared  and  presented  in  an  international  meeting  and  national 
meetings . 

In  the  past,  we  have  analyzed  flow  over  missiles  with  cru- 
ciform fins  in  subsonic  and  supersonic  Mach  Number  ranges  (Ref- 
erences 1,  2,  3).  The  above  includes  studies  on  the  effect  of 
viscosity  and  viscous  interaction  of  the  flow  on  the  jet  plume. 
All  the  results  reported  hereto  were  for  clean  configurations. 
During  the  present  contract  period  an  effort  was  made  to  study 
the  flow  over  configurations  with  "steps"  or  "cutouts".  This 
was  done  using  slenderbody  theory.  The  results  of  this  study 
are  reported  in  detail  in  Section  II  of  this  report. 

We  have  studied  the  fin-body  combinations  at  lower  angles 
of  attack  up  to  10°  in  the  past.  Due  to  the  growing  interest 
in  higher  angles  of  incidence  flight,  efforts  were  made  to 
study  fin  body  combinations  at  moderately  high  angles  of  attack. 
For  fins,  the  wake  has  to  be  modelled  effectively.  This  is  done 
by  letting  the  vortex  filaments  in  the  wake  to  float,  and  to 
decide  their  own  course.  The  forebody  at  these  angles  of  attack 
generates  two  vortices  which  are  symmetrical.  The  effects  of 
these  vortices,  (the  strengths  of  which  are  given  by  slenderbody 
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theory),  are  also  included  in  the  computations.  The  detailed 
investigations  and  results  are  reported  in  Section  III  of  this 
report . 

There  has  been  a growing  interest  exhibited  in  high  angle 
of  attack  flows.  The  vortices  generated  by  the  forebody  become 
unsymmetrical  at  this  point,  and  due  to  the  shedding  of  vortices, 
the  large  and  small  vortices  switch  sides.  This  results  in 
switching  of  directions  of  side  forces  which  are  of  the  order  of 
the  normal  force.  These  effects  were  studied  by  extending  our 
method  for  the  ogive  and  conical  forebodies.  These  results  are 
reported  in  Section  IV. 

During  the  present  contract  period  three  papers  were  pub- 
lished. The  titles  are  presented  below,  and  a full  copy  of  the 
papers  is  included  in  appendices  to  this  report.  The  three  papers 
are : 

1.  "An  Analysis  of  Various  Fin  Geometries  and  Their 
Influence  on  Rocket  Performance",  by  N.  Uchiyama, 

T.  H.  Moulden  and  J.  M.  Wu,  presented  at  Proceedings 
of  the  Twelfth  International  Symposium  on  Space 
Technology  and  Science,  Tokyo,  1977,  pp.  109-120. 

2.  "An  Analysis  of  Wing-Body  Combinations  at  Moderate 
Angles  of  Attack",  (AIAA  Paper  Number  78-62),  by 
N.  Uchiyama,  R.  P.  Mikkilineni,  and  J.  M.  Wu,  at 
AIAA  16th  Aerospace  Sciences  Meeting,  Huntsville, 

Alabama,  1978. 

3.  "Study  of  Asymmetric  Nose  Vortices  of  Bodies  of 
Revolution",  (AIAA  Paper  Number  79-27),  by  Y.  Niwa, 

R.  P.  Mikkilineni,  and  J.  M.  Wu,  at  17th  Aerospace 
Sciences  Meeting,  New  Orleans,  Louisiana,  1979. 
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II.  FLOW  PAST  A CUT  OUT 

Flow  over  bodies  with  axisymmetric  cut  outs  was  investigated, 
and  the  results  are  reported  in  this  section.  In  the  case  of 
axisymmetric  bodies  with  a clean  configuration  (i.e.,  no  cut  outs 
etc.,)  the  cylindrical  portion  of  the  body  has  an  almost  flat 
pressure  distribution.  Looking  at  some  experimental  results  (Ref- 
erences 4 and  5)  in  Figures  1 and  2,  we  clearly  see  the  effect 
the  cavity  has  on  the  nearly  sonic  flow.  The  results  reported  in 
those  figures  are  for  two-dimensional  flows.  If  the  cut  out  is 
very  shallow,  it  is  nearly  a two-dimensional  flow.  As  results 
for  three-dimensional  or  axisymmetric  flows  were  not  accessible 
at  the  time  of  this  investigation,  two-dimensional  results  were 
used  in  studying  the  nature  of  the  flow. 

Looking  at  Figures  1 and  2 (and  some  more  experimental  res- 
ults) , it  can  be  seen  that  there  are  re-circulatory  flow  regions 
near  the  two  edges,  longer  near  the  front  edge.  It  is  hence  pos- 
sible to  consider  a stream  line  that  is  closest  to  the  body  and 
passes  over  the  re-circulatory  flow  regions  (treating  those  reg- 
ions more  as  a part  of  the  solid  body) . (Figures  3 and  4) . This 
technique,  it  can  be  seen,  would  avoid  the  detailed  viscous  flow 
calculations  and  the  shock-boundary  layer  interaction  studies. 

The  main  approximation  made  in  the  analysis  is  that  the  flow 
enters  the  cut  out  at  an  angle  of  about  7°  and  leaves  the  cut  out 
at  about  45° . The  results  from  these  calculations  are  reported 
in  subsequent  figures  (Figures  5 - 9) . It  can  be  seen  that  the 
nature  of  the  pressure  distribution  agrees  pretty  well  with  the 
experimental  results. 

Also,  the  pressure  distribution  over  the  cavities  was  found 
to  be  a function  of  the  three  non-dimensional  parameters,  namely, 
Mach  Number,  Reynold's  Number,  and  the  height  to  length  ratio  of 
the  cavity  (h/L) . It  is  possible  to  compute  the  streamline  shape 
(Figure  5)  from  the  experimental  pressure  distribution,  and  an 
equation  can  be  formulated  to  the  approximated  streamline  in  terms 
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of  the  three  non-dimensional  parameters  mentioned  above.  This 
technique  is  essentially  computing  the  equivalent  body  in  the 
flow.  The  advantage  of  this  method  is  that  once  the  equations 
are  formulated,  it  is  possible  to  compute  the  pressure  distri- 
bution at  other  values  of  the  non-dimensional  parameters  using 
the  equivalent  body  approach  mentioned  in  the  first  part  of 
the  section.  Such  a study  was  undertaken  during  the  present 
contract  period,  but  it  was  not  completed. 
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III.  FIN-BODY  COMBINATION  AT  MODERATELY 
HIGH  ANGLE  OF  ATTACK 

A computational  technique  is  developed  to  study  the  aero- 
dynamic flow  over  wing  (or  fin)  body  combinations  at  moderately 
high  angles  of  attack.  Numerical  results  including  the  static 
stability  derivatives  and  the  damping  in  pitch  and  roll  are  pre- 
sented for  a typical  missile  configuration  with  fins  and  can- 
ards . 

The  flow  over  a fin-body  combination  at  moderate  angles  of 
attack  is  sketched  in  Figure  10.  One  of  the  main  characteris- 
tics of  this  flow  is  the  separated  flow  at  the  fin  tips  and/or 
the  leading  edge.  The  vortex  sheet  separated  from  these  sharp 
edges  rolls  into  strong  vortices  above  the  wing.  The  strong 
cross  flow  produced  by  these  vortices  is  responsible  for  larger 
lift  coefficients  of  the  low  aspect-ratio  wings  at  high  angles 
of  attack  (Reference  6).  Further,  this  interaction  introduces 
non-linearities  into  the  problem.  For  angles  of  attack  up  to 
about  25  degrees,  the  body  supports  two  symmetric  vortices 
originating  from  the  nose.  For  these  angles  of  attack,  the  flow 
over  the  body  is  not  separated  on  the  leeward  side,  and  conse- 
quently, the  asymmetric  vortex  system  of  high  angle  of  attack 
flow  over  slender  bodies  need  not  be  considered. 

The  bound  vortex  sheet  of  the  fins  is  represented  by  a vor- 
tex lattice.  This  is  done  by  dividing  the  fin  into  a finite 
number  of  elements  each  carrying  a bound  vortex  along  the  span 
and  two  trailing  vortices  along  the  chord.  The  separated  sheets 
from  the  edges  and  tips  are  modelled  by  a series  of  non-inter- 
secting vertical  filaments.  Each  filament  in  turn  is  represen- 
ted by  a series  of  straight  line  vortex  segments  joined  together 
at  the  points  called  nodes.  This  representation  is  shown  in 
Figure  11. 

The  body  surface  is  divided  into  a finite  number  of  rec- 
tangular elements  with  each  element  carrying  a concentrated 
singularity  of  source-sink  type  at  its  centroid.  This  represen- 
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cation  assumes  that  the  flow  is  not  separated  on  the  leeward 


side  of  the  body.  This  is  shown  in  Figure  12.  The  representa- 
tion of  the  nose  is  not  straight  forward.  One  has  to  take  into 
consideration  the  potential  vortices  generated  by  the  lifting 
nose  in  addition  to  the  surface  distribution  of  source  type 
singularity.  In  this  work,  the  surface  of  the  nose  is  repre- 
sented by  using  concentrated  source-sink  type  singularities  at 
the  centroids  of  the  rectangular  elements  representing  the  sur- 
face. The  potential  vortices  are  represented  based  either  on 
the  theoretical  results  or  empirical  data.  (Reference  7) . The 
vortices  are  represented  by  a horse  shoe  vortex  system,  the 
strength  and  location  of  which  are  predetermined  from  empirical 
relations  (Reference  7) . The  representation  of  the  nose  is 
shown  in  Figure  13. 

The  total  velocity  on  the  body  is  the  resultant  of  the 
contributions  from: 

a)  free  stream 

b)  body  rotation  in  pitch  and  roll 

c)  singularity  distribution  on  the  body;  and 

d)  the  free  vortices 

The  boundary  condition  on  the  body  to  be  satisfied  is  that  the 
normal  component  of  this  resultant  velocity  be  zero.  Further, 
the  free  vortices  in  the  vortex  wake  are  force  free  i.e.,  they 
are  aligned  with  the  local  velocity  vector.  The  unknowns  in 
this  representation  are: 

a)  the  strength  of  the  concentrated  singularities 

b)  the  strength  of  the  bound  vortices  on  the  wing 
surface,  and 

c)  the  orientation  of  the  free  vertical  filaments  in 
the  stream. 

These  unknowns  are  to  be  determined  by  the  boundary  conditions 
on  the  body  and  the  condition  on  the  free  vertical  elements . 

As  this  problem  is  non-linear,  the  unknowns  are  determined 
using  an  iterative  scheme  which  includes: 
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a)  Initial  direction  and  strength  are  assigned  to  each 
vortex  element  in  the  wake. 

b)  With  the  wake  fixed  the  strength  of  the  vortices  and 
sources  on  the  body  are  determined  by  requiring  that 
the  resulting  flow  satisfies  the  boundary  condition 
on  each  element  of  the  body  and  the  wing. 

c)  With  the  strength  of  the  singularities  fixed  the 
wake  is  re-oriented  to  satisfy  the  force-free  condi- 
tion on  the  vertical  filaments.  Steps  (b)  and  (c) 
are  repeated  until  the  iterates  converge. 

Using  the  technique  discussed  in  earlier  sections,  different 
configurations  of  typical  missiles  are  analyzed.  Results  are  pre- 
sented in  graphical  form  for  two  configurations.  The  first  con- 
figuration shown  in  Figure  14  is  analyzed  for  the  angle  of  attack 
a = 19.4  degrees. 

As  the  flow  is  symmetric  about  x-z  plane,  this  property  can 
be  used  to  reduce  the  number  of  computations.  The  representation 
is  only  needed  for  half  of  the  configuration.  The  discretization 
involves  the  division  of  the  fin  into  5x5  elements,  the  body 
into  6 x 12  panels. 

The  convergence  of  the  iteration  scheme  is  shown  in  Figure 
15. 

The  normal  force  coefficient  is  plotted  as  a function  of  the 
iteration  number.  It  is  seen  that  the  process  converged  within  8 
iterations . 

In  Figure  16,  the  spanwise  load  distribution  is  shown.  The 
vortex  system  in  the  plane  perpendicular  to  the  axis  of  the  body 
and  containing  the  trailing  edge  is  shown  in  Figure  17. 

The  rolled  up  vortex  from  the  wing  is  clearly  identified  in 
this  figure.  Further,  the  nose  vortex  is  shown  in  this  figure. 

The  nose  vortex  is  a concentrated  vortex  located  above  the  body 
about  three  diameters  away.  The  roll-up  of  the  separated  vortices 
is  shown  in  Figure  18.  In  Figure  18a,  the  formation  is  shown  at 
half  chord.  Figure  18b  shows  the  vortex  at  the  trailing  edge,  and 
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Figure  18c  shows  the  vortex  in  the  wake,  one  chord  length  away 
from  the  trailing  edge. 

Figures  19  and  20  show  the  trailing  vortex  filaments  in  x-y 
and  x-z  planes.  These  figures  clearly  show  the  rolling  up  of  the 
vortical  filaments. 

The  second  configuration  studied  using  this  technique  is 
shown  in  Figure  21.  Compared  to  the  first  configuration,  this 
has  higher  aspect  ratio  wing.  The  flow  is  analyzed  for  an  angle 
of  attack  of  19.4  degrees. 

The  convergence  of  the  normal  force  coefficient  is  shown  in 
Figure  22.  The  process  converged  within  8 iterations.  Figure 
21  shows  the  variation  of  the  load  along  the  span  of  the  wing. 

In  Figure  24  the  formation  of  the  vortex  due  to  separation 
on  the  wing  is  shown  at  stations  x = 0.5,  x = 1.0,  and  x = 2. 

In  the  wake  corresponding  to  x = 2 , a clearly  formed  vortex  can 
be  seen.  Figure  16  shows  the  vortex  system  in  a plane  perpendi- 
cular to  the  axis  of  the  missile  at  x = 1.0. 

The  modelling  and  computation  of  a flow  over  fin-body  com- 
binations for  any  angle  of  attack  and  orientation  is  not  an  easy 
task.  In  this  work,  a computational  technique  is  described  and 
used  to  compute  the  flow  over  a practical  missile  configuration 
at  moderately  high  angles  of  attack.  Although  the  results  are 
presented  for  rectangular  fins  only,  the  scheme  could  effectively 
be  used  for  any  other  planform.  At  these  angles  of  attack,  the 
linear  theories  are  not  suitable  because  of  the  nonlinear  behavior 
of  the  flow. 

For  a detailed  mathematical  formulation  as  well  as  computa- 
tional procedures  and  results,  see  Appendix  2 of  this  report. 

Further  knowledge  and  work  are  needed,  especially  on  the 
viscous  flow  behavior  including  flow  separation,  to  effectively 
treat  the  fin-body  combinations  at  higher  angles  of  attack.  This 
is  due  to  the  asymmetric  flow  over  the  slender  bodies  at  angles 
of  attack  larger  than  25°.  Any  reasonable  computational  model 
should  incorporate  the  aspect  of  the  asymmetric  flow  in  the  over- 
all representation. 
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IV.  NOSE  VORTICES  AT  HIGH  INCIDENCE 


I 


The  existence  of  asymmetric  vortex  flow  on  a body  of  revo- 
lution at  high  angles  of  attack  or  at  incidence  with  zero  side 
slip  to  any  plane  due  to  the  symmetricity  of  the  body  has  been 
observed  a long  time  ago  (Reference  8) . Until  recently,  this 
asymmetric  vortex  flow  phenomenon  was  not  paid  much  attention 
because  flying  at  very  high  angle  of  incidence  is  rarely  done. 

This  situation  has  been  changed,  however.  The  increasing  demand 
for  better  performance  of  missiles  has  necessitated  new  improve- 
ments in  characteristics  of  missiles  at  high  angles  of  attack, 
which  can  occur  during  an  extremely  sharp  maneuver.  The  asym- 
metric vortices  generate  an  induced  side  force  and  yawing  mom- 
ent when  the  body  is  at  high  angles  of  attack  and  has  zero  side 
slip.  The  side  force  and  yawing  moment  may  force  the  missile 
into  either  an  uncontrollable  or  an  unpredictable  situation. 

The  worse  case  results  in  the  sudden  switching  of  direction  of 
the  side  force.  Much  of  the  research  work  on  the  study  of  asym- 
metric vortex  flow,  both  theoretical  and  experimental,  has  been 
done  because  improvement  of  aerodynamic  characteristics  in  this 
area  has  been  required. 

Excellent  background  reviews  on  high  angle  of  attack  tech- 
nology have  been  given  in  References  9-11.  One  of  the  basic 
approaches  of  using  cross  flow  analogy  by  appealing  to  a two- 
dimensional  analysis  was  originally  proposed  by  Bryson  (Reference 
12)  and  Angelucci  (Reference  13),  which  is  followed  by  this  study. 
The  method  used  here  is  the  Wardlaw's  modified  version  (Reference 

14)  of  cross-flow  analogy  to  an  asymmetric  system.  It  is  inter- 
esting to  note  here  that  Sarpkaya's  impulsive  analogy  (Reference 

15) ,  Marshall  and  Def fenbaugh ' s inclusion  of  boundary-layer  con- 
tribution (Reference  16) , are  also  helpful  to  our  understanding  of 
lee- side  vortex  flow. 

The  main  object  of  this  study  is  a first-cut  estimation  on 
the  behavior  of  two  concentrated  asymmetric  vortices  formed  near 
the  forebody.  Fidler,  et.  al  (Reference  17)  indicated  that  the 


I 
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vortex  flow  may  not  have  a concentrated  vortex  core  but  a diffused 
and  distributed  vorticity.  In  this  first-cut  study,  however,  it 
is  assumed  that  the  vortex  core  exists  in  the  front  portion  of  the 
forebody.  Many  flow  field  visualization  experiments  support  this 
notion  (See,  for  instance,  Reference  18).  An  attempt  is  then  made 
to  compute  the  normal  force  and  the  side  force  magnitudes  as  well 
as  traces  of  this  pair  of  asymmetric  vortices.  The  previous  paper 
of  Uchiyama,  Mikkilineni  and  Wu  (Reference  19)  has  shown  that  the 
first  pair  of  vortices  generated  by  any  body  of  revolution  with 
pointed  nose  is  the  key  to  leading  the  subsequent  vortices  from 
the  main  body.  Their  work,  however,  has  been  addressed  to  the 
symmetric  pair  of  vortices  and  applicable  only  up  to  a moderately 
high  incidence  case. 

A symmetric  vortex  flow  is  generally  observed  on  the  lee- 
ward side  of  bodies  of  revolution  in  the  subsonic  to  moderate 
supersonic  flow  range  at  moderate  angles  of  attack  and  asymmetric 
vortex  flow  is  observed  at  high  angles  of  attack.  The  schematic 
flow  fields  on  these  two  cases  are  shown  in  Figure  25.  The  vor- 
tex systems  in  a cross  flow  plane  are  shown  in  Figures  26-27. 

The  two  conditions  used  in  solving  the  motion  completely 
are  as  follows:  The  first  of  the  conditions  is  that  the  feeding 
sheet  and  the  roll-up  concentrated  vortex  must  have  zero  net 
force  on  the  free  vortex  system.  So,  the  concentrated  vortex 
itself  must  experience  an  equal  but  opposite  force  as  generated 
by  the  feeding  sheet.  The  second  is  similar  to  the  Kutta  condi- 
tion in  that  the  location  of  the  stagnation  points  where  the  flow 
leaves  the  body  has  to  be  assigned.  They  are  taken  here  as  points 
where  the  wall  sheer  stress  vanishes. 

Governing  equations  as  derived  previously  are  integrated 
numerically  by  using  the  Runge-Kutta  routine  with  step  size  con- 
trol . 

The  integration  of  wake  vortex  flow  requires  starting  from 
the  boundary  layer  flow  separation  point.  However,  it  is  very 
difficult  to  set  up  this  point.  Bryson  (Reference  12)  and 
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Angelucci  (Reference  13)  used  0q  = 40  degrees  measured  from  the 
horizontal  plane  as  shown  in  Figure  3 for  the  separation  angle. 
Wardlaw  used  42  degrees  and  Clark  and  Peoples  (Reference  20) 
used  65  degrees.  In  this  study  40  degrees  has  been  chosen  for 
separation  angle.  A study  has  been  conducted  to  investigate  the 
effect  of  this  separation  angle.  The  influence  of  this  angle  on 
the  normal  force  was  found  to  be  small,  but  a large  change  in  the 
side  force  was  observed. 

A symmetric  vortex  system  which  has  the  strength  of  Agt  is 
calculated  first.  After  that  an  asymmetrioity  by  using  a free 
parameter  is  instroduced.  This  a symme tricity , used  in  this  study, 
is  given  by  one  of  the  following  two  methods.  In  the  first  method, 
a vortex  position  is  artificially  moved  toward  the  body  and  another 
vortex  is  moved  outward  from  the  body  in  half  of  the  asymmetricity , 
respectively.  This  procedure  has  been  designated  as  "Case  1"  in 
the  subsequent  calculation.  In  the  second  method,  it  has  been 
arbitrarily  assumed  that  one  vortex  keeps  its  same  position  while 
the  other  vortex  is  moved  outward  from  the  body  in  an  asymmetric 
pattern.  This  case  has  been  designated  "Case  2"  in  the  present 
calculation.  Case  1 and  Case  2 are  illustrated  in  Figure  28. 

Three  values  of  asymmetricities  of  0.5  percent,  1 percent  and  2 
percent  of  local  radius  are  studied.  The  ogive  nose  and  the  cone 
nose,  both  with  a fineness  ratio  of  3.5,  are  chosen  in  this  study. 
The  flow  chart  of  numerical  computation  procedure  is  given  in 
Figure  29. 

The  computer  results  for  the  ogive  forebody  and  the  cone- 
nose  cases  are  discussed  below.  The  results  for  the  two  cases 
are  in  general  similar;  however,  some  distinct  differences  also 
exist . 

1.  Ogive  Nose 

The  computed  normal  force  component  acting  on  an  ogive  fore- 
body of  3.5  fineness  ratio  and  at  various  angles  of  attack  are 
shown  in  Figures  30  and  31.  The  Figure  30  computation  is  based 


on  the  asymmetricity  starting  condition  Case  1 as  described  pre- 
viously. The  Figure  31  computation  is  based  on  the  starting  con- 
dition of  Case  2.  The  two  results  are  very  similar  over  the 
entire  range  of  the  angles  of  attack  studied.  There  is  a very 
minute  difference  in  the  normal  force  coefficient  values.  Case 
2 predicts  a slightly  larger  normal  force  than  Case  1.  The  cor- 
responding circulation  for  Case  2 is  also  slightly  larger  than 
Case  1.  For  the  Case  1,  the  magnitude  of  total  circulation 
(l*i I + I *2 I)  is  very  close  to  the  symmetrical  vortex  flow  case 
by  using  a modified  slender  body  theory  as  proposed  by  Fidler  and 
Bateman  (Reference  21) . The  circulation  generated  by  the  Case  2 
procedure  is  slightly  larger  than  the  symmetrical  vortex  flow 
case.  For  Case  1 and  angle  of  attack  below  30  degrees  computa- 
tion, it  was  observed  that  the  smaller,  or  closer-to-the-body 
inner  vortex  has  a tendency  first  to  move  in  and  then  move  out 
in  the  radial  direction.  For  Case  2,  the  two  vortices  have  been 
shown  to  move  out  monotonically  in  all  of  the  computations. 

As  far  as  in  the  reality,  it  is  not  known  how  the  asymmetri- 
city is  introduced  into  the  flow  field.  The  values  of  free  para- 
J meter  used  in  indicating  the  percentage  of  starting  asymmetricity 

are  shown  in  the  above-mentioned  two  figures.  At  this  stage  we 
j have  no  way  of  knowing  the  correct  value  of  asymmetricity  corres- 

ponding to  the  real  situation.  In  Figure  30,  two  experimental 
results  by  Coe,  et  al.  (Reference  18)  and  Chapman  and  Keener 
(Reference  22)  are  also  included  for  comparison.  It  is  inter- 
esting to  note  that  these  two  sets  of  data  do  not  agree  with  each 
other  in  spite  of  the  data  having  been  taken  over  similar  config- 
urations and  with  the  same  model  size.  Our  computed  results  are 
closer  to  Chapman's  data  than  Coe's  measurement,  especially  when 
the  angles  of  attack  are  below  45  degrees.  At  larger  angles  of 
attack  the  present  calculated  results  depart  from  Chapman's  data 
and  approach  Coe ' s measurement . 

The  relationship  between  the  computed  side  force  and  angles 
of  attack  is  shown  in  Figure  32  for  Case  1 and  Figure  33  for 
Case  2.  How  it  is  observed  that  the  Case  1 procedure  generates 
a slightly  higher  side  force  than  the  Case  2 procedure  in  this 
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computation.  This  is  just  opposite  to  the  normal  force  case. 
However,  the  difference  is  still  very  minute.  The  magnitude  of 
side  force  is  more  susceptible  to  the  asymmetricity  than  the 
normal  force.  The  curves  are  spread  more  as  the  angle  of  attack 
increases.  Both  Chapman's  and  Coe's  data  are  plotted  for  compar- 
ison with  the  computed  results.  At  lower  angles  of  attack  (say 
below  45  degrees),  Coe's  measurement  lies  in  the  low  asymmetri- 
city curves  while  Chapman's  data  shows  a higher  asymmetricity. 

Our  computation  showed  that  in  order  to  achieve  such  a high  side 
force,  a twice-larger  asymmetric  vortices  strength  has  to  be 
assumed  to  start  the  computation  (see  Figure  34) . 

The  computed  results  on  the  outer  vortex  flow  along  the 
x-axis  for  different  angles  of  attack  are  shown  in  Figure  35. 

The  initial  asymmetricity  in  the  radial  distance  is  given  as  two 
percent  (Case  2)  in  this  calculation.  It  is  found  that,  up  to 
a = 35  degrees,  the  strength  of  this  vortex  flow  remains  near  a 
constant  (with  a slight  increase)  as  the  vortex  flow  moves  toward 
downstream.  As  the  angle  of  attack  exceeds  35  degrees,  the  vortex 
strength  increases  as  it  moves  toward  downstream  direction.  This 
indicates  more  feeding  of  the  vorticity  to  the  vortex  sheet  occurs 
at  a higher  angle  of  incidence.  This  is  because  the  cross-flow 
component  gains  magnitude  as  the  U-sina  increases.  However,  the 
increase  in  the  strength  is  in  a nonlinear  way  which  indicates 
the  contribution  of  the  dr/dx  term  is  highly  coupled  with  the  a 
influence.  From  this  study,  it  is  also  realized  that  the  growth 
in  asymmetricity  is  nonlinear  for  very  high  incidence  cases. 

A comparison  of  the  inner  versus  the  outer  vortex 
strength  has  been  made.  The  computed  non-dimensional  vortex 
strengths  of  A^  and  A2  along  the  x-axis  are  shown  in  Figure  36. 

The  asymmetricity  was  taken  as  two  percent  and  a = 55  degrees  in 
this  calculation.  The  growth  rate  of  these  two  vortices  is  quite 
different  initially,  but  at  downstream  it  becomes  nearly  the  same. 
It  is  found  that  the  strength  of  these  two  vortices  has  a tendency 
to  equalize  in  their  magnitudes,  but  the  asymmetricity  in  their 
distance  becomes  larger  and  larger.  The  strength  ratios  of  these 
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two  vortices  along  the  x-axis  are  shown  in  Figure  37.  As  the 
incidence  angle  a is  less  than  35  degrees,  the  strength  ratio 
remains  nearly  the  same;  however,  it  increases  rapidly  for  a>40 
degrees.  Note  the  difference  in  the  growth  rate  along  the  x-axis 
for  the  different  ot  in  Figure  27.  Also,  note  the  decreasing 
trend  of  the  strength  ratio  as  a increases.  Finally,  the  vortex 
strength  ratio  approaches  unity  as  a limit  for  very  high  a case. 

Typical  calculated  traces  for  the  inner  and  the  outer  vor- 
tex cores  along  the  x-axis  are  given  in  Figure  38.  For  this  com- 
putation, it  was  assumed  that  a = 45  degrees,  and  the  asymmetri- 
city  of  radial  distance  is  0.5  percent.  The  result  indicated  that 
two  vortices  become  farther  apart  as  they  flow  toward  downstream. 

2.  Conical  Forebody 

A similar  computation  has  been  carried  out  for  the  case  of 
a conical  forebody.  This  forebody  has  the  same  fineness  ratio 
of  3.5  as  mentioned  previously. 

The  calculated  normal  force  components  at  various  angles  of 
incidence  are  shown  in  Figure  39.  For  the  conical  forebody  case, 
it  appears  that  the  influence  of  asymmetricity  is  much  smaller 
than  in  the  ogive  nose  case.  A comparison  with  Coe's  data  is  also 
included  in  Figure  39.  The  theory  predicts  rather  well  against 
this  set  of  data  up  to  a = 45  degrees.  At  higher  a,  the  discrep- 
ancy shows  up.  It  is  difficult  to  explain  this  discrepancy  at 
higher  incidence  by  only  considering  the  variation  in  normal 
force  component. 

The  computed  side  force  components  against  the  variation  in 
angles  of  attack  are  shown  in  Figure  40.  The  agreement  of  data 
with  the  theory  is  rather  poor.  Especially,  the  measurement 
indicated  a switching  in  the  side  force  direction  at  a = 40  degrees, 
which  cannot  be  predicted  by  the  present  calculation.  In  the  pre- 
vious ogive  case,  it  has  been  shown  a reasonable  agreement  with 
limited  sets  of  data  on  the  side  force  components.  This  indicates 
that  the  side  force  component  prediction  is  much  more  difficult 
for  the  conical  forebody.  A drastic  change  in  the  vortex  pattern 
must  occur  at  ot  = 40  degrees.  The  prediction  of  switching  of  dir- 
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ection  in  the  side  force  component  remains  to  be  a formidable 
problem  at  this  time. 

The  computed  strength  of  outer  vortex  along  the  x-axis  is 
shown  in  Figure  41.  The  strength  of  the  vortex  is  increased 
very  quickly. 

The  calculated  strength  of  inner  and  outer  vortices  along 
the  x-axis  is  shown  in  Figure  42.  In  the  conical  forebody  case, 
both  vortices  are  developed  very  quickly  compared  to  that  of  the 
ogive  case.  It  is  interesting  to  observe  that  the  strength  of 
the  inner  vortex  exceeds  the  outer  vortex  as  the  flow  moves  tow- 
are  downstream.  Again,  a good  experiment  is  needed  to  confirm  or 
deny  such  a phenomenon. 

A similar  condition  as  shown  in  Figure  38  for  the  ogive  case 
(namely  a = 45  degrees  and  0.5  percent  asymmetricity)  has  been 
applied  to  compute  traces  of  outer  and  inner  vortices  for  the 
cone  case.  This  result  is  shown  in  Figure  43.  A marked  differ- 
ence in  traces  is  observed  between  these  two  nose  configurations, 
especially  in  the  outer  vortex.  In  the  case  of  the  conical  nose, 
the  outer  vortex  moves  out  significantly  away  from  the  body  com- 
pared to  that  of  the  ogive  case  --  as  shown  in  Figure  38.  The 
distance  between  the  inner  and  the  outer  vortices  is  larger  in 
the  conical  forebody  case  as  well  (see  Figure  43) . For  both  the 
cone  and  the  ogive  cases,  the  inner  vortex  core  remains  very 
close  to  the  body.  It  is  interesting  to  observe  that  the  strengths 
of  the  two  vortices  (Figure  42)  remain  about  the  same  in  spite  of 
large  differences  in  the  distance  between  them.  In  the  real  flow 
case,  it  is  probably  very  difficult  to  maintain  a very  extended 
vortex  feeding  sheet  from  the  body  surface  to  the  core.  A break- 
down of  this  feeding  sheet  results  in  a vortex  shedding.  Conse- 
quently, a third  vortex  may  appear  near  the  body  surface.  There 
is  good  reason  to  suspect  that  this  third  and  close-to-the-body 
vortex  is  the  cause  of  side  force  switching. 

Finally,  a comparison  of  traces  or  positions  of  vortices 
for  the  ogive  and  the  cone  is  made  and  is  given  in  Figure  44. 

The  positions  shown  in  Figure  44  are  viewed  from  the  front  of  the 
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body.  A significant  difference  between  each  of  the  vortex  traces 
is  illustrated.  In  the  ogive  case,  the  trace  is  somewhat  curved, 
while  in  the  cone  case,  the  trace  is  nearly  a straight  line.  A 
noticeable  difference  between  the  ogive  and  cone  is  also  clear. 

It  may  conclude  that  the  forebody  radius  distribution,  i.e.,  dr/dx, 
affects  rather  strongly  the  motion  of  vortex  flow.  In  spite  of 
the  small  asymmetrical  value  in  radius  initially  assigned,  the 
final  result  in  asymmetry  is  quite  large,  especially  for  the  cone 
case.  A check  against  Foppl  theory  (Reference  23)  is  made.  Foppl 
line  is  the  locus  of  centers  of  two  symmetric  vortices  on  the  lee- 
side  of  a two-dimensional  cylinder.  Our  result  showed  that  the 
path  of  vortices  does  not  follow  the  Foppl  line,  which  was  also 
observed  by  Fidler,  et  al.  (Reference  17). 

The  detailed  formulation  of  computational  scheme  is  given  in 
Appendix  2. 


V.  CONCLUDING  REMARKS 


The  above  remarks  summarize  the  work  performed  under  the 
present  study.  As  a result  of  this  work,  it  is  evident  that 
several  fundamental  questions  still  need  to  be  resolved. 

1.  Potential  flow  analysis  is  adequate  for  predicting 
forces  if  the  strength  and  location  of  the  free- 
vortices  is  known. 

2.  The  initial  formation  of  the  vortex  on  the  body 
(either  a forebody  or  a lifting  planar  surface) 
is  a viscous  flow  problem. 

3.  Work  is  needed  to  study  this  viscous/ inviscid 
flow  interaction. 
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FIGURE  1.  Experimental  Result  On  Shallow  Cavity  (Ref.  5) 
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FIGURE  8.  Plot  of  Section  Lift  Coefficient  Along  Body  Surface 
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FIGURE  9.  Plot  of  Pressure  Coefficient  Along  Body  Surface 
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FIGURE  9.  (Continued) 
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FIGURE  10.  Schematic  flow  field  around  missile  configuration 

at  moderately  high  angle  of  attack 
with  emphasis  on  fin  effects. 
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FIGURE  15.  Variation  of  CN  with  iteration 

number  (configuration). 
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FIGURE  22. 
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Variation  of  Cn  with  iteration 
number  (configuration  II). 
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Figure  27.  Asymmetric  Vortex  System  in 
a Cross  Flow  Plane. 


Figure  30 . Normal  Force  Coefficient 
Attack  for  Ogive  Forebody 
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Normal  Force  Coefficient  vs.  Angles 
Attack  for  Ogive  Forebody  (Case  2) . 


Strength  of  Outer  Vortex  Along  x-Axis 
Ogive  Forebody  (Case  2) . 


Strength  of  Inner  and  Outer  Vortex 
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Figure  37.  The  Vortex  Strength  Ratio  Along  x-Axis  (Case  1) . 


Figure  38.  The  Trace  of  the  Inner  and  Outer 
Vortex  Along  x-Axis  (Case  2) . 
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Figure  39.  Normal  Force  Coefficient  vs.  Angies  of 
Attack  for  Conical  Forebody  (Case  2). 
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Figure  41.  The  Absolute  Strength  of  Outer 
Vortex  Along  x-Axis  (Case  2) . 


Along  x-Axis  (Case  2) . 


